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Abstract 

An  improved  one-dimensional  dynamic  model  of  a  tubular  SOFC  stack  capable  of  system  simulation  in  the  virtual  test  bed  (YTB)  simulation 
environment  is  presented  in  this  paper.  This  model  is  based  on  the  electrochemical  and  thermal  modeling,  accounting  for  the  voltage  losses  and 
temperature  dynamics.  The  modeling  of  an  external  reformer  is  also  included  in  this  study.  A  detailed  parametric  analysis  of  working  conditions 
and  cell  configuration  of  the  solid  oxide  fuel  cell  (SOFC)  stack  is  the  main  focus  of  this  paper.  The  following  operating  parameters  are  investigated: 
pressure  ratio,  temperature,  mass  flow  rate,  external  reforming  degree  and  stream  to  carbon  (S/C)  ratio.  The  cell  geometric  parameters  studied 
include  cell  diameter  and  cell  length.  Elevated  operating  pressure  improves  the  cell  performance.  Whereas,  higher  operating  temperature  decreases 
both  the  Nernst  potential  and  the  irreversible  losses,  resulting  in  an  initial  increase  then  a  decrease  in  cell  efficiency.  It  was  found  that  a  higher  S/C 
ratio  yields  a  lower  H2  concentration  and  partial  pressure,  which  has  a  negative  effect  on  the  Nernst  potential.  Increased  cell  diameter  is  found  to 
increase  the  power  due  to  a  larger  activation  area  at  the  same  time  and  due  to  longer  current  path  length  there  is  an  increase  in  the  ohmic  loss. 
Increased  length  of  the  cell  has  the  undesired  affect  of  an  increased  pressure  drop. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  high  temperature  solid  oxide  fuel  cell  (SOFC),  with  the 
advantages  of  its  high  efficiency,  environmental  friendliness  and 
flexibility  of  usable  fuel  types,  has  been  considered  as  one  of 
the  most  promising  technologies  for  electric  energy  generation. 
Many  variations  of  SOFC  fuel  cell  design  are  possible:  such  as 
tubular,  planar  and  monolithic.  The  tubular  SOFC  power  system 
(Siemens- Westinghouse  Power  Corporation),  integrated  with  a 
gas  turbine,  has  achieved  an  overall  55%  efficiency  [1].  The 
planar  and  monolithic  SOFC  technologies  are  frequently  ana¬ 
lyzed  in  the  literature.  The  current  paper  focuses  on  an  improved 
dynamic  model  for  a  tubular  SOFC  stack,  and  a  detailed  param¬ 
eter  analysis  of  operating  conditions  and  cell  geometric  config¬ 
urations.  Parameter  analysis  not  only  shows  the  influence  of  the 
parameters’  variance  on  the  SOFC  performance,  but  also  iden¬ 
tifies  areas  of  the  largest  possible  impact  with  the  least  effort. 
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In  the  past  couple  of  years,  several  models  have  been  devel¬ 
oped  and  tested  to  study  the  design  and  operating  conditions 
of  SOFC  stack.  Most  models,  however,  are  steady  state  mod¬ 
els  [2-4]  and  are  valid  for  only  the  specific  operating  points. 
Padulles  et  al.  [5]  developed  a  SOFC  model  with  focus  on 
electrochemistry  and  species  dynamics,  but  temperature  and 
heat  transfer  dynamics  were  not  considered  in  their  model.  A 
dynamic  transient  SOFC  model  has  been  developed  by  Sedghisi- 
garchi  [6],  however,  their  model  is  based  on  lumped  capacitance 
model  where  the  spatial  distributions  of  temperature,  current 
density  and  species  concentration  are  not  investigated.  In  addi¬ 
tion,  few  studies  have  evaluated  the  effect  of  length  of  current 
pathway,  which  has  significant  influence  on  the  cell  diameter. 

In  this  paper,  a  one-dimensional  dynamic  model  of  a  tubular 
SOFC  is  presented  with  external  reforming  and  system  inte¬ 
gration  capabilities.  Based  on  the  electrical  behavior,  chemical 
reaction  equilibrium  and  energy  balance,  the  proposed  model 
can  predict  the  SOFC  characteristics  in  steady  states  and  also 
in  transient  operating  states.  The  accuracy  and  reliability  of  the 
model  is  demonstrated  by  comparisons  with  experimental  data 
from  the  literature.  The  distinctive  feature  of  the  current  model 
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Nomenclature 

cp  specific  heat  (J  (mol-1  K-1)) 

Dh  hydraulic  diameter  (cm) 

E  activation  energy  (82  kJ  mol- 1 ) 

E°  Nemst  potential  (V) 

F  Faraday  constant  (96487  C  mol- 1 ) 

G  Gibbs  free  energy  (J  mol- 1 ) 

h  enthalpy  (J  mol-1) 

1  current  density  (A  cm-2) 

z'o  exchange  current  density  (A  cm-2) 

iL  limiting  current  (A  cm-2) 

k  pre-exponential  factor  (4274) 

^shift  equilibrium  constant  for  the  shifting  reaction 

p  partial  pressure  (bar) 

m  molar  number  of  control  volume  (mol) 

hi  molar  flow  rate  of  species  i  (mol  s-1) 

Nu  Nusselt  expression 

^ch4  conversion  rate  for  CH4  per  unit  area  of  the  cell 
(mol(s-1  cm-2)) 

R  Universal  gas  constant  (8.314  J  (mol-1  K-1)) 

t  time  (s) 

T  absolute  temperature  (K) 

Q  heat  transfer  (kJ) 

Uf  fuel  utilization 

x  conversion  rates  of  CH4  (mol  s-1) 

y  conversion  rates  of  CO  (mol  s-1) 

2  conversion  rates  of  FU  (mol  s  - 1 ) 

Greek  symbols 

a  electron  transfer  coefficient 

0  efficiency 

77 cone  concentration  loss  (V) 

rj Act  activation  loss  (V) 

Ohm  ohmic  resistance  loss  (V) 

p  electrical  resistance  (£2) 

Subscripts 
a  air 

ar  air  in  reaction 

ap  preheated  air 

i  inlet 

e  outlet 

f  fuel 

s  solid  part 


is  that  the  model  takes  into  consideration  the  variation  of  vari¬ 
ables  in  the  axial  direction,  as  a  result  it  can  predict  the  cell 
performance  more  accurately.  Another  important  feature  of  the 
proposed  model  is  that,  the  concentration  loss  is  modeled  by 
introducing  partial  pressure  related  to  limiting  current  and  an 
equivalent  circuit  has  been  built  to  evaluate  how  the  length  of 
the  current  pathway  affects  the  ohmic  loss  of  a  tubular  SOFC.  In 
addition,  the  simulation  results  of  the  model,  which  is  dynamic 
in  nature,  can  be  used  for  performance  evaluation  and  cell  design 


optimization  under  variable  operating  conditions  and  geometric 
conditions. 

The  virtual  test  bed  (VTB)  platform  provides  an  effective 
computational  environment  to  simulate  the  dynamic  perfor¬ 
mance  of  the  SOFC  stack  [7,8].  The  non-linear  model  equations 
based  on  electrochemistry  and  thermodynamics  are  discretized 
in  resistive  companion  (RC)  form  for  effective  implementation 
in  the  VTB  platform. 

The  paper  is  organized  as  follows.  In  Section  2,  the  system 
configurations,  operating  conditions  and  important  assumptions 
of  SOFC  stack  are  presented.  The  model  description  is  presented 
in  details  in  Section  3.  The  parameter  analysis  of  operating 
conditions  and  cell  configuration  are  presented  in  Section  4. 
Conclusions  are  drawn  in  Section  5. 

2.  System  configuration 

Fig.  1  illustrates  the  system  configurations  of  a  SOFC  stack 
implemented  in  the  VTB  simulation  environment.  The  system 
consists  of  one  SOFC  stack,  an  external  reformer,  two  valves, 
one  electrical  load  and  two  thermal  sinks.  A  methane-steam 
mixture  is  supplied  to  the  external  reformer  before  being  deliv¬ 
ered  to  the  SOFC  stack.  The  inlet  pressure  and  temperature  of 
fuel  and  air  are  assumed  as  constant  and  defined  as  user  input 
parameters.  The  fuel  from  the  external  reformer  and  the  air 
are  channeled  to  the  anode  and  cathode  of  SOFC,  respectively. 
Hydrogen  diffuses  through  the  anode  and  reacts  with  oxygen 
ions.  In  the  cathode,  oxygen  diffuses  through  the  electrode  and 
the  free  electrons  are  absorbed.  Electrical  energy  is  hence  pro¬ 
duced  along  with  the  generation  of  the  heat.  The  electrical  energy 
is  supplied  to  the  electrical  load.  Under  the  assumption  of  adi¬ 
abatic  conditions,  a  portion  of  the  heat  generated  during  the 
SOFC  electrochemical  reaction  and  SOFC  internal  resistance  is 
consumed  to  reform  the  fuel,  and  the  rest  is  used  to  heat  up  the 
gases.  For  simplicity,  the  exhaust  gases  from  SOFC  stack  are 
channeled  to  the  thermal  sinks  directly  in  the  system  modeled 
here.  It  can  be  noted  that  this  high  temperature  gas  is  a  prime  tar¬ 
get  for  energy  recovery,  i.e.,  the  high  temperature  exhaust  gases, 
which  consist  of  unutilized  fuel  and  depleted  air,  can  be  supplied 
to  a  combustor  for  operation  of  a  Brayton  cycle  gas  turbine.  The 
detailed  operating  conditions  and  SOFC  parameters  used  in  the 
current  study  are  listed  in  Table  1 . 

With  a  given  inlet  pressure,  the  mass  flow  rates  of  the  fuel 
and  air  are  controlled  by  adjusting  the  openings  of  two  throt¬ 
tle  valves.  In  the  external  reformer,  methane  participates  in  the 
reforming  reaction,  which  produces  hydrogen  and  carbon  diox¬ 
ide.  The  electrical  load  is  simplified  as  a  resistive  load,  whose 
value  is  treated  as  an  input  parameter  and  can  be  changed  to 
study  the  behavior  of  the  SOFC  under  variable  load  conditions. 
The  SOFC  model  shown  in  Fig.  1  is  a  two-level  VTB  model.  The 
first  level,  shown  as  a  SOFC  stack  icon,  includes  one  electrical 
type  connector  and  four  fluid  type  connectors,  which  are  an  air 
inlet,  an  air  outlet,  a  fuel  inlet  and  a  fuel  outlet  respectively.  It 
can  be  noted  that  each  of  the  fluid  type  connector  consists  of 
two  ports:  one  natural  port  and  one  signal  port.  The  natural  port 
is  used  for  the  variable  pressure  and  mass  flow  rate  while  the 
signal  port  is  used  for  temperature  and  species  fraction.  More 
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Table  1 

Operating  conditions  and  cell  specifications 


Parameter 

Value 

Unit 

Operation  conditions 

Stack  power  output 

440 

kW 

Cell  voltage 

0.603 

V 

Average  current  density 

0.301 

Acm~2 

Average  power  density 

0.182 

W  cm-2 

Air  pressure  ratio 

5 

NU 

Fuel  pressure  rato 

5 

NU 

Air  flow  rate 

1.08 

kgs-1 

Fuel  flow  rate 

0.118 

mols-1 

Air  inlet  temperature 

900 

K 

Fuel  inlet  temperature 

900 

K 

Stack  pressure  loss 

2 

% 

Fuel  utilization 

0.855 

NU 

Fuel  composition 

33.3%  CH4 

66.7%  H20 

Cell  parameters 

Number  of  cell 

10000 

Cell  area 

242 

cm  x  cm 

Cell  length 

50 

cm 

electrolyte  thickness 

0.004 

cm 

Anode  thickness 

0.19 

cm 

Cathode  thickness 

0.01 

cm 

information  on  port  definition  can  be  found  in  reference  [7].  By 
connecting  the  ports  of  SOFC  with  those  of  the  models  available 
in  VTB  model  library,  user  can  easily  integrate  the  SOFC  model 
into  any  specific  system.  Fig.  2  illustrates  the  second  level  of  the 
SOFC  model,  which  includes  the  detailed  configurations  of  the 
singular  cell.  As  shown  in  Fig.  2,  the  cell  is  divided  into  several 
elements  along  the  flow  direction  by  the  finite  volume  method 
in  which  all  the  governing  equations,  i.e.  electrochemical  reac¬ 
tion,  mass  conservation  and  energy  conservation  equations,  are 
solved  for  each  element.  The  number  of  elements  is  determined 
by  the  complexity  and  accuracy  requirement  of  the  system  and 
can  be  easily  changed  by  the  user  as  needed.  The  detailed 


description  of  the  finite  volume  method  is  presented  in  Section 
3.3. 

3.  Model  description 

A  one-dimensional  dynamic  model  of  a  SOFC  stack,  based 
on  the  earlier  work  [20],  has  been  developed.  The  stack  model 
has  been  improved  significantly  by  introduction  of  an  external 
reformer,  appropriate  modeling  of  concentration  loss  by  intro¬ 
ducing  partial  pressure  related  to  limiting  current,  and  most 
importantly,  an  equivalent  circuit  has  been  built  to  evaluate  how 
the  length  of  current  pathway  affects  the  ohmic  loss  of  tubu¬ 
lar  SOFC.  The  following  important  assumptions  are  made  in 
developing  the  SOFC  model. 

1.  One-dimensional  behavior  along  the  stream  direction. 

2.  Walls  of  the  SOFC  tube  and  external  walls  of  the  SOFC 
system  enclosure  are  assumed  adiabatic.  There  is  no  axial 
heat  exchange  between  cells. 

3 .  Gas  distribution  among  cells  is  uniform.  There  is  no  deviation 
of  the  gas  distribution  among  cells. 

4.  No  gas  phase  reactions  occur. 

5.  No  gases  leak  outside  of  the  system. 

6.  All  gases  are  assumed  to  be  ideal  gas. 

3.1.  External  reformer  model 

Two  reactions  occur  in  the  external  reformer. 

Reforming  CH4  +  F^CXA)  CO  +  3H2 

(Ah,  =  206.1  kJ  mol-1)  (1) 

Shifting  CO  +  H20(  v)  +»  C02  +  H2 

(Ahs  =  —41.1  kJ  mol-1)  (2) 


Fig.  2.  Second  level  of  the  SOFC  model. 
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Table  2 

Equilibrium  constant  for  Eq.  (5) 


A 

B 

C 

D 

E 

5.47E-12 

— 2.57E— 08 

4.64E-05 

— 3.93E— 02 

1.32E+01 

where  x  and  y  are  the  conversion  rates  of  CH4  and  CO,  respec¬ 
tively.  The  conversion  rate  x  is  calculated  from  inlet  molar  flow 
rate  of  methane  and  degree  of  external  reforming. 

X  =  degr  •  «ch4  (3) 


where  hlc H4  is  reformer  inlet  molar  flow  rate  of  methane  and 
degr  is  degree  of  external  reforming,  which  is  given  as  an  input 
parameter. 

The  shifting  reaction  is  considered  to  reach  thermodynamic 
equilibrium,  which  can  be  calculated  from  the  partial  pressures 
of  the  reactants  and  products. 


^shiftcn 


^H2  •  ”co2 
^h2o  '  ^co 


The  equilibrium  constant  for  the  shifting  processes  is  tem¬ 
perature  dependent  and  can  be  expressed  as  a  polynomial  Eq. 

(5) 


log  ATghift  =  AT4 +  BT3  +  CT2  +  DT  +  E  (5) 


where  the  values  of  the  coefficients  [10]  are  listed  in  Table  2. 
The  mass  balance  for  the  external  reformer  is  given  as: 


"ch4  =  "ch4  —  X 

(6) 

«h2  =  «h2  +3  x  +  y 

(7) 

”h2o  =  "h2o  -  x  —  y 

(8) 

«co  =  «co  +  x  ~  y 

(9) 

«co2  =  «co2  +  >’ 

(10) 

The  energy  balance  for  the  external  reformer  is  given  to  be: 


y  hihiTi )  =  y  heh(Te)  +  x  •  A h,  +  y  ■  A hsb  (1 1) 


3.2.  SOFC  electrochemical  model 


Fig.  3  illustrates  the  configuration  and  current  pathway  of 
a  single  tubular  cell.  Hydrogen  diffuses  through  the  anode  and 
reacts  with  oxygen  ions.  In  the  cathode,  oxygen  diffuses  through 
the  electrode  and  the  free  electrons  are  absorbed.  The  electrolyte 
can  conduct  oxygen  ions  at  temperatures  above  800  °C.  The 
anode  of  one  cell  is  connected  to  the  cathode  of  the  next  via  the 
interconnect.  Physical  dimensions  of  each  of  the  components 
are  shown  in  Table  1 . 

The  concerned  reactions  inside  the  cell  can  be  summarized 
as  follows 

Reforming  CH4  +  H20(a)  CO  +  3H2 

(Ah,  =  206.1  kJ  mol-1)  (12) 


Tubular  SOFC 


Interconnect 


Fuel  Electrode 
Electrode 

Air  Electrode 

Air  Feed  Tube 


Electron  flow 
(current) 


Fig.  3.  Configuration  and  flow  path  of  electrons  and  ions  in  a  cross-section  of  a 
tubular  SOFC  tube. 


Shifting  CO  +  H20(y)  C02  +  H2 

(Ahs  =  -41.1  kJ  mol"1)  (13) 


Electrochemical  H2  +  yCMi)  — »■  H2O 

(A/iec  =  -241.8  kJ  mol-1)  (14) 


where  x ,  y  and  z  are  the  conversion  rates  of  CH4,  CO  and 
H2,  respectively.  A  hY,  A  hs  and  A  hQC  stand  for  reaction  heat  of 
reforming,  shifting  and  electrochemical  reactions,  respectively. 
The  reforming  reaction  is  a  highly  endothermic  reaction.  The 
conversion  rate  x  is  assumed  to  be  controlled  by  a  conversion 
rate  per  unit  area  of  the  cell,  which  is  named  fcH4  given  by 
Achenbach  [9]. 


^ch4  =  k  •  pcn4  •  exp 


—E 


R  •  T 


(15) 


where  k  is  the  pre-exponential  factor  with  a  value  of  4274,  E  the 
activation  energy  of  82  kJ/mol,  and  pcw4  is  the  partial  pressure 
of  methane. 

Similar  to  the  external  reforming  process,  under  the  assump¬ 
tion  that  shifting  reaction  is  always  in  equilibrium,  the  equilib¬ 
rium  constants  can  be  calculated  from  Eq.  (4). 

Hydrogen  conversion  rate  is  directly  related  to  the  current, 
by  Faraday’s  law. 

z=I--F  (16) 

The  fuel  utilization  coefficient  Uf  is  given  by 


Uf  = 


z 


hlU2+h 


co 


+  4  hi 


ch4 


(17) 


The  cell  electrical  power  is  calculated  as  the  product  of  cell 
current  and  the  cell  electrodes  voltage.  The  reversible  voltage 
calculated  based  on  the  knowledge  of  the  cell  potential  from  the 
Nernst  equation. 


—  A  G° 
IF 


+ 


RT 

2F 


In 


1/2 

Ph2  Pq2 


P  h2o 


(18) 


where  p  terms  are  the  partial  pressures  of  reacting  species,  R  the 
universal  gas  constant  and  F  is  Faraday  constant.  The  change  of 
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Gibbs  free  energy  A  G°  is  the  function  of  temperature  and  can 
be  described  as  a  polynomial  equation.  The  coefficients  a-c  are 
taken  from  Reid  et  al.  [11] 

A  G°=aT2+bT  +  C  (19) 

The  Nernst  potential  is  reduced  due  to  the  following  irre¬ 
versible  losses: 


1.  concentration  loss  Oconc  (V); 

2.  activation  loss  77 Act  (V); 

3.  ohmic  resistance  loss  ijohm  (V). 


Therefore,  the  actual  voltage  can  be  expressed  as 
V  =  E  (l) Act  T  OOhm  T  OConc )  (20) 


(1)  The  concentration  loss  is  calculated  as  follows  [12,13]: 


RT  (  i\ 

OConc  —  —  In  (  1  —  T"  ) 

nF  \  iLJ 


(21) 


where  i  is  the  cell  current,  /l  the  limiting  current  and  n  is  the 
number  of  electrons  participating  in  the  reaction.  Account¬ 
ing  for  the  known  impact  of  pressure  on  mass  diffusion  and 
experimental  observation,  variable  z'l  can  be  written  as  [19]: 


il  =  ilx 


(22) 


where  a  is  a  constant,  p  the  pressure  and  “primed”  variables 
are  for  atmospheric  conditions. 

(2)  The  activation  loss  occurs  when  the  rate  of  electrochemical 
reaction  at  an  electrode  surface  is  controlled  by  sluggish 
electrode  kinetics  [14,15],  which  is  shown  as  follows 


RT  f  Gen 
7?Act  =  — -  In  — 
ant  \  iq 


(23) 


where  a  is  the  electron  transfer  coefficient  of  the  reaction  at 
the  electrode  and  io  is  the  exchange  current  density,  which 
is  a  function  of  stack  temperature. 

(3)  According  to  Bessett  [16],  the  resistivity  for  the  different 
layers  of  tubular  SOFC  is  given  by 


p  =  0.008 114  e600/ T  (air  electrode) 
p  =  0.00294  e10350//r  (electrolyte) 
p  =  0.00298  e_1392//r  (fuel  electrode) 
p  =  0 . 1 25 68  e4690/ T  (interconnect) 


(24) 

(25) 

(26) 
(27) 


The  current  pathway  in  a  tubular  SOFC  is  shown  in  Fig.  3. 
The  current  collects  circumferentially  around  the  anode  and  the 
cathode  tubes;  and  oxide  ions  move  perpendicularly  across  the 
electrolyte  layer.  To  account  for  the  ohmic  loss  along  the  cir¬ 
cumferential  pathway,  an  equivalent  circuit  shown  in  Fig.  4  was 
created  to  evaluate  the  ohmic  loss.  The  ohmic  loss  for  the  elec¬ 
trolyte  layer  and  interconnect  can  be  given  by 

R  =  P  ■  j  (28) 


Fig.  4.  Equivalent  circuit  for  the  tubular  SOFC  cross. 


where  L  is  the  layer  thickness  and  A  is  unit  area.  But  for  the  anode 
and  cathode  layers,  the  effect  of  circumferential  current  pathway 
must  be  considered  to  evaluate  the  ohmic  loss.  The  length  of  cur¬ 
rent  pathway  is  the  function  of  cell  diameter.  The  influence  of 
cell  diameter  is  investigated  by  simulating  the  equivalent  circuit 
in  Fig.  4  with  different  diameters.  Therefore  a  corrective  coef¬ 
ficient  was  introduced  to  account  for  anode  and  cathode  layer 
ohmic  losses. 


(29) 


where  cl  is  a  constant,  D  the  cell  diameter  and  D'  is  the  reference 
diameter. 

The  ohmic  loss  can  be  expressed  as 


OOhm  —  I  -  R 


(30) 


3.3.  SOFC  thermal  model 

Fig.  5  shows  the  schematic  configuration  of  a  single  tubular 
fuel  cell  with  discretized  axial  elements  used  for  heat  transfer 
calculations.  The  heat  released  from  chemical  reaction  and  elec¬ 
trical  resistance  loss  is  absorbed  by  the  solid,  where  it  is  further 
transferred  to  the  fuel  and  air  steam.  The  air  entering  from  the 
air  feed  tube  is  pre-heated  by  the  air  in  reaction.  For  the  thermal 
model  the  following  assumptions  are  made 

1.  Heat  release  and  absorption  arising  from  reforming,  shift¬ 
ing,  electrochemical  reactions  and  electrical  resistance  occur 
within  solid  part  of  each  cell. 

2.  Heat  transfer  between  the  solid  and  the  gas  steams  occur 
by  convection  heat  transfer.  Radiation  heat  transfer  between 
solid  part  and  gas  steams  is  not  considered  here.  Preliminary 
calculation  by  the  authors  revealed  that  radiation  heat  transfer 
was  much  smaller  than  the  convective  heat  transfer,  as  a  result 
radiation  heat  transfer  was  neglected.  The  effect  of  radiation 
will  be  investigated  in  more  detains  in  the  future  work. 

3.  The  axial  conduction  heat  transfer  between  nodes  is 
neglected. 

4.  The  cell  voltage  at  each  element  is  uniform. 
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Fig.  5.  Heat  transfer  assumption  and  element  dividing. 


As  seen  in  Fig.  5,  the  cell  is  divided  into  elements  along  the 
flow  direction.  The  temperature  gradient  of  gas  steams  in  the 
flow  direction  is  assumed  to  depend  only  on  the  convection  heat 
transfer  between  the  channel  walls  and  the  gas  stream.  For  each 
element,  as  shown  in  Fig.  6,  there  are  four  control  volumes:  fuel, 
solid,  reaction  air  and  pre-heated  air.  Four  different  temperatures 
7f,  Ts,  Tar,  Tap,  corresponding  to  the  temperatures  of  fuel,  solid, 
air  in  reaction  and  pre-heated  air,  respectively,  are  considered 
within  any  element  of  each  cell  under  the  assumption  of  uniform 
temperature  within  each  of  the  control  volume. 

For  the  solid  control  volume,  the  energy  equation  is  given  as 

Ms  •  Cp  •  =  (2 gen  -Qf-Qar-W  (31) 

where  Qgen  is  the  total  internal  heat  generation,  Gf  and  Gar  the 
convection  heat  transfers  to  fuel  and  reaction  air,  respectively 
and  W  is  the  electrical  power  produced  by  the  SOFC  stack. 

The  total  heat  generation  Ggen  accounting  for  the  reactions 
(shifting,  reforming  and  electrochemical)  and  resistance  to  cur¬ 
rent  is  given  by 

(2 gen  =  —  (x  •  A hT  +  y  •  A/zsh  +  z  •  A hec)  +  72  •  ^  Ri  (32) 

The  electrical  power  is  calculated  as  the  product  of  current 
and  voltage 


W  =  V  •  /  •  ACeii 


(33) 


Ts 

Tar 

Tap 


The  convection  heat  transfers  Gf  and  Gar  can  be  calculated 
by 

Gf  =  ocfAf(Ts  -  Tf )  (34) 

Gar  —  °Gr  Aar{Ts  —  Tar)  (35) 


The  heat  transfer  coefficients  o?f  and  aar  are  given  by  the 
Nusselt  expression: 


Nu  •  kc 

a  =  - 

Dh 


(36) 


where  Nusselt  number  is  set  to  4  according  to  Rokni  and  Yuan 
[17],  kc  is  the  conductivity  of  the  gas  and  Dh  is  hydraulic  diam¬ 
eter. 

The  validity  of  numerical  scheme  used  for  calculating 
the  instantaneous  temperature  in  the  gas  control  volume  is 
based  on  a  small  simulation  time  step.  For  a  small  simula¬ 
tion  time  step  (0.0001  s),  the  originally  simultaneous  process 
of  the  flow  of  the  stream  flowing  and  the  transfer  of  heat  can 
be  assumed  to  be  two  separate  processes,  one  preceding  the 
other: 


(1)  Firstly,  stream  flows  along  the  axial  direction,  as  seen  in 
Fig.  7,  at  a  given  time  step  some  gas  flows  into  the  control 
volume  while  an  amount  of  gas  flows  out  into  the  next  ele¬ 
ment.  We  assume  each  of  the  control  volume  to  be  well 
mixed  (this  is  a  valid  assumption  for  small  control  vol¬ 
ume),  we  can  now  calculate  an  average  temperature  using  an 
energy  balance.  This  process  can  be  described  by  Fig.  7. The 
resulting  conservation  of  energy  equation  is  given  by 

h  •  At 

hj  =  hr  T  \hr  —  hr  ]  •  -  (37) 

^ave  %  1  L  -ZonJ  v  / 

m 

where  7]n,  Tor\  and  Tave  are  the  temperatures  of  entering 
steam,  exiting  steam  and  the  well  mixed  control  volume, 


Fig.  6.  Control  volume  defining  for  one  element. 


Fig.  7.  Gas  control  volume. 
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Fig.  8.  Comparison  of  experimental  data  [18]  with  model,  parametric  pressure 
variation  fuel:  89%  H2,  11%  H2O,  85%  fuel  utilization;  operating  temperature: 
1000 °C. 


•  Temperature,  including  stack  mean  temperature,  air  inlet  tem¬ 
perature  and  fuel  inlet  temperature. 

•  Mass  flow  rate,  including  air  flow  rate  and  fuel  flow  rate. 

•  Degree  of  external  reforming. 

•  Stream  to  carbon  (S/C)  ratio. 

The  cell  geometric  parameters  studied  include: 

•  diameter  of  the  tubular  SOFC,  and 

•  length  of  the  tubular  SOFC. 

It  should  be  pointed  out  that  all  the  calculations  are  made  at 
the  design  point  with  following  input  parameters:  the  tempera¬ 
ture  of  fuel  and  air,  the  flow  rate  of  fuel  and  air  and  cell  geometry, 
which  are  kept  constant  unless  mentioned  specifically. 


respectively,  h  the  molar  flow  rate,  At  is  the  simulation  time 
step  and  m  is  the  molar  number  of  control  volume. 

(2)  Heat  transfer  occurs  after  the  mixing.  The  energy  equations 
for  gaseous  control  volumes  can  be  calculated  by 


fuel  mfCp 


dTf 
d  t 


dTar 

reaction  air  marcD - 

p  d  t 


(38) 

(39) 


preheated  air  rafCp 


dTap 
d  t 


(40) 


The  SOFC  stack  efficiency 
the  electric  power.  LHV  is 
fuel. 


is  given  by  Eq.  (41),  where  W  is 
the  lower  heating  value  of  inlet 


^Overall  — 


w 

LHV 


4.  Parameter  analysis 


(41) 


The  parameter  analysis  of  tubular  SOFC,  based  on  the  design 
point  listed  in  Table  1  and  system  shown  in  Fig.  1,  will  be  pre¬ 
sented  in  this  section.  The  following  operating  parameters  are 
investigated: 


•  Stack  pressure  ratio  (operation  pressure  to  atmosphere  pres¬ 
sure). 


4.1.  Pressure  ratio 

Fig.  8  shows  the  comparisons  between  simulation  results  and 
experimental  data  [18]  for  I-V  characteristics  curve  at  the  oper¬ 
ating  pressures  ratio  of  1,  5  and  10.  A  very  good  agreement  is 
obtained  throughout  the  test  range.  For  the  sake  of  comparing 
with  the  experimental  conditions,  the  fuel  composition  is  set  to 
89%  H2,  11%  H2O,  with  85%  fuel  utilization.  The  mean  oper¬ 
ating  temperature  is  set  at  1000  °C  by  adjusting  the  air  flow  rate. 

Under  the  simulation  conditions  listed  in  Table  2,  the  effects 
of  operating  pressure  on  the  cell  voltage  and  efficiency  is  shown 
in  Fig.  9(a).  As  operating  pressure  increases,  the  cell  voltage 
and  efficiency  increase  because  of  the  increased  Nemst  poten¬ 
tial  and  lower  concentration  loss.  When  operating  pressure  ratio 
increases  from  1  to  15,  an  improvement  of  7.38%  and  24.17%  is 
obtained  for  the  cell  voltage  and  efficiency,  respectively.  How¬ 
ever,  the  improvement  of  cell  voltage  and  efficiency  is  not  a  linear 
function  of  the  increase  of  pressure.  After  an  initial  increase  the 
increase  is  less  with  further  increase  in  pressure.  This  is  due  to 
the  fact  that  the  pressure  contribution  to  the  Nernst  potential  is 
logarithmic  in  nature.  The  solid  temperature  profile  under  vari¬ 
ous  operating  pressures  is  shown  in  Fig.  9(b).  It  is  obvious  that 
the  temperature  profile  becomes  sharper  with  higher  pressure. 
The  temperatures  of  the  solid  near  exit  are  nearly  the  same  for 
different  operating  pressure  while  there  is  a  big  increase  in  max¬ 
imum  temperature  for  the  higher  pressure.  Given  a  constant  fuel 
and  air  mass  flow  rate,  a  higher  pressure  leads  to  a  better  cell 
efficiency  and  fuel  utilization,  resulting  in  a  bigger  amount  of 
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Fig.  9.  Effect  of  stack  pressure  on  the:  (a)  cell  voltage  and  efficiency;  (b)  solid  part  temperature  profile. 
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Current  Density  (A/cm2) 

Fig.  1 0.  Comparison  of  experimental  data  [2 1  ]  with  model,  parametric  operating 
temperature  variation  fuel:  89%  H2,  11%  H2O,  with  85%  fuel  utilization. 

fuel  taking  part  in  the  reaction  in  the  inlet.  As  a  result,  there  is  a 
steeper  rise  in  the  temperature  for  the  higher  pressure  due  to  the 
fact  that  more  heat  is  generated 

4.2.  Temperature 
4.2.1.  Stack  temperature 

The  mean  stack  temperature  is  adjusted  and  kept  to  1173, 
1223  and  1273  K  by  changing  the  air  mass  flow  rates  in  order 
to  validate  the  comparison  with  the  experimental  data  [21].  The 
fuel  composition  for  this  comparison  is  89%  H2, 11%  H2O,  with 
85%  fuel  utilization. 

As  shown  in  Fig.  10,  the  agreement  is  very  good  (<1.5%)  for 
the  whole  1173  K  case  and  good  for  the  lower  current  density 


part  of  the  1223  K  case.  For  the  1273  K  case  and  exit  side  for 
the  1223  K  case,  relatively  higher  variance  may  be  due  to  an 
overestimation  of  the  activation  loss  and  the  simplification  of 
the  SOFC  geometry.  However,  the  trend  of  the  simulated  I-V 
characteristics  matches  well  with  that  of  the  experimental  data. 
Moreover,  the  maximum  relative  discrepancy  between  the  model 
and  the  experiment  results  are  below  5%,  showing  a  fairly  good 
fit.  Along  with  the  comparison  for  pressure  in  Section  4. 1 ,  it  can 
be  validated  that  the  proposed  SOFC  model  indeed  works  well 
to  predict  the  overall  performance  of  a  tubular  SOFC. 

4.2.2.  Air  inlet  temperature 

The  effect  of  air  inlet  temperature  on  the  SOFC  performance 
is  shown  in  Fig.  11  (a  and  b).  The  air  inlet  temperature  is  varied 
and  all  the  other  input  parameters  are  kept  constant.  As  a  matter 
of  fact,  the  impact  of  air  inlet  temperature  variance  on  SOFC  per¬ 
formance  is  caused  by  the  solid  part  temperatures  changes  due  to 
the  varied  air  inlet  temperature.  As  shown  in  Fig.  1 1  (b),  the  mean 
solid  temperature  increases  from  1180  to  1390K  when  the  air 
inlet  temperature  rises  from  600  to  800  K.  The  reduction  of  irre¬ 
versible  losses  while  operating  under  higher  temperature  leads  to 
a  better  performance  in  terms  of  cell  voltage,  Uf  and  efficiency. 
Therefore  the  voltage,  Uf  and  efficiency  rise  at  first.  However, 
this  effect  is  completely  cancelled  out  by  a  further  reduction  of 
Nernst  potential  because  of  the  higher  temperature.  Therefore, 
cell  voltage,  Uf  and  efficiency  increase  initially  as  the  air  temper¬ 
ature  is  increased  and  then  decrease  with  the  increase  of  air  inlet 
temperature  (when  the  decrease  in  Nernst  potential  is  greater 
than  the  corresponding  reduction  in  the  irreversible  losses). 


Fig.  11.  Effect  of  air  inlet  temperature  on  the:  (a)  cell  voltage  and  Uf,  (b)  mean  temperature  and  efficiency. 


(a)  Fuel  inlet  Temperature  (k) 


(b)  Fuel  inlet  Temperature  (k) 


Fig.  12.  Effect  of  fuel  inlet  temperature  on  the:  (a)  cell  voltage  and  Uf,  (b)  mean  temperature  and  Efficiency. 
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Fig.  13.  Effect  of  air  flow  rate  ratio  on  the:  (a)  cell  voltage  and  Uf  (b)  mean  temperature  and  efficiency. 


Fig.  14.  Effect  of  fuel  flow  rate  ratio  on  the:  (a)  cell  voltage  and  Uf  (b)  mean  temperature  and  efficiency. 


4.2.3.  Fuel  inlet  temperature 

The  effect  of  fuel  inlet  temperature  on  the  SOFC  performance 
is  shown  in  Fig.  12(a  and  b).  Because  of  a  relatively  lower  mass 
flow  rate  of  fuel,  there  is  a  minor  influence  on  the  cell  perfor¬ 
mance  when  the  fuel  inlet  temperature  is  changed.  As  shown  in 
Fig.  12(b),  there  is  only  a  total  increase  of  68  K  in  the  mean  tem¬ 
perature  when  fuel  inlet  temperature  rises  from  600  to  800  K.  In 
this  small  range  of  operating  temperature,  the  effect  of  reduction 
in  irreversible  losses  overcomes  the  influence  of  the  reduction  in 
Nemst  potential.  So  the  cell  performance  in  terms  of  cell  volt¬ 
age,  Uf,  and  efficiency  is  improved  with  the  increasing  of  fuel 
inlet  temperature. 

4.3.  Mass  flow  rate 
4.3.1.  Airflow  rate 

Fig.  13(a  and  b)  show  the  curves  of  cell  voltage,  Uf,  mean 
solid  temperature  and  efficiency  via  the  air  flow  rate  ratio.  The 


simulation  begins  from  the  operating  condition  listed  in  Table  2. 
The  other  input  parameters  are  kept  constant  while  air  flow  rate 
is  increased  by  10%,  20%,  50%  by  adjusting  the  opening  of 
the  control  valve  in  the  testing  system.  The  solid  temperature 
decreases  because  of  the  cooling  effect  of  larger  amount  of  air. 
The  decreased  temperature  yields  to  a  lower  cell  voltage,  Uf  and 
lower  cell  efficiency. 

4.3.2.  Fuel  flow  rate 

Similar  to  the  effect  of  increasing  air  flow  rate,  higher  fuel 
flow  rate  yields  a  decreased  mean  solid  temperature  as  shown  in 
Fig.  14(b).  With  increasing  of  fuel  flow  rate,  Uf,  and  efficiency 
decreases  quickly  as  shown  in  Fig.  14(a  and  b);  this  is  because 
fuel  can  only  be  partially  oxidized.  The  cell  voltage  increases 
initially  because  of  a  higher  amount  of  fuel  participating  in  the 
reaction.  But  this  effect  will  be  overcome  later  due  to  the  cooling 
effect  of  higher  fuel  flow  rate.  Therefore,  cell  voltage  increase 
at  first  and  then  decrease  with  the  increase  of  fuel  flow  rate. 


Fig.  15.  Effect  of  external  reforming  on  the  SOFC  performance,  (a)  Cell  voltage  and  Uf,  (b)  mean  temperature  and  efficiency. 
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Fig.  16.  Effect  of  S/C  ratio  on  the  cell  voltage  and  Uf. 


Fig.  18.  Effect  of  length  of  tubular  SOFC  on  the  cell  voltage. 


4.4.  Degree  of  external  reforming 

The  effect  of  degree  of  external  reforming  on  the  SOFC  per¬ 
formance  is  shown  in  Fig.  15(a  and  b);  with  a  degree  variation 
from  10%  to  90%.  The  mass  flow  rate  of  fuel  and  air  and  the 
inlet  temperature  of  fuel  and  air  are  kept  constant.  The  effect  of 
internal  reforming  is  weak  with  increasing  the  degree  of  exter¬ 
nal  reforming  because  less  and  less  methane  will  be  reformed 
inside  of  the  SOFC.  With  increasing  degree  of  external  reform¬ 
ing,  there  is  less  CH4  undergoing  reforming  inside  of  SOFC  and 
accordingly  less  endothermic  reforming  reactions.  As  a  result, 
the  mean  solid  temperature  increases.  The  cell  performance  in 
terms  of  cell  voltage,  Uf,  and  efficiency  benefits  from  the  higher 
mean  solid  temperature,  as  shown  in  Fig.  15 (a  and  b). 

4.5.  S/C  ratio 

The  effect  of  stream  to  carbon  (S/C)  ratio  on  the  SOFC  per¬ 
formance  is  shown  in  Fig.  16.  The  S/C  ratio  is  adjusted  from 
1.5  to  4,  but  the  mass  flow  rate  of  methane  is  kept  constant  by 
adjusting  the  opening  of  the  control  valve.  As  shown  in  Fig.  16, 
the  cell  voltage  and  Uf  decreases  due  to  the  increase  of  S/C 
ratio.  The  reason  for  this  lies  in  the  fact  that  the  higher  S/C  ratio 
yields  a  lower  H2  concentration  and  partial  pressure,  which  has 
negative  effect  on  the  Nernst  potential. 

4.6.  Diameter  of  the  tubular  SOFC 

Three  different  diameters:  D=15.88,18  and  22  mm  are  inves¬ 
tigated.  The  diameter  of  the  tubular  SOFC  is  adjusted  by  chang¬ 


ing  the  diameter  of  air  inlet  tube  listed  in  Table  2.  The  thicknesses 
for  different  layers  are  kept  constant.  Uf  and  mean  solid  temper¬ 
ature  are  fixed  throughout  the  analysis.  As  suggested  earlier  in 
Section  3.2  regarding  ohmic  losses,  increasing  the  diameter  of 
tubular  SOFC  leads  to  an  undesired  increase  in  the  current  path 
length  which  in  turn  yields  a  higher  ohmic  loss,  resulting  in  an 
decrease  in  the  Nernst  potential  of  the  cell.  The  diameter  of  the 
SOFC  has  less  effect  on  the  activation  and  concentration  losses. 
Fig.  17(a)  shows  the  effect  of  diameter  variation  on  the  cell  volt¬ 
age.  Obviously,  for  a  given  current  density,  a  higher  cell  voltage 
can  be  obtained  with  a  smaller  diameter.  As  shown  in  Fig.  17(b), 
a  larger  diameter  SOFC  produces  a  higher  output  power  due  to 
the  increased  cell  area.  However,  the  power  density  of  a  larger 
diameter  case  is  smaller  than  that  of  smaller  diameter  case. 

4. 7.  Length  of  tubular  SOFC 

The  effect  of  cell  length  on  the  SOFC  performance  is  shown 
in  Fig.  18.  Three  cases  are  investigated,  with  L  =  50,  100  and 
150  cm,  respectively.  Similar  to  the  diameter  analysis,  the  thick¬ 
nesses  for  different  layers  are  kept  constant.  Uf  and  mean  solid 
temperature  are  fixed  throughout  the  analysis.  The  three  V-I 
curves  overlap  with  one  another,  which  reveals  that,  unlike  the 
variation  in  diameters,  length  has  no  effect  on  the  cell  voltage 
for  a  given  current  density  because  it  does  not  change  the  cur¬ 
rent  path  length.  But  the  total  power  output  increases  because 
of  the  larger  cell  reaction  area.  In  the  simulation,  the  pressure 
drop  inside  the  SOFC  is  given  and  kept  constant.  However,  in 
reality,  the  effect  of  pressure  drop  inside  the  SOFC  is  enhanced 
due  to  the  higher  mass  flow  rate.  So  the  diameter  of  cell  has 


D=1 5.88  mm  — • —  D=1 8  mm  — ■ —  D=22  mm 


Fig.  17.  Effect  of  diameter  of  tubular  SOFC  on  the:  (a)  cell  voltage;  (b)  power. 
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to  be  increased  to  counteract  the  pressure  drop  effects,  which 
will  have  the  undesired  effect  of  increased  current  path  length 
resulting  in  increased  cell  ohmic  resistance. 

5.  Conclusions 

A  one-dimensional  dynamic  SOFC  model,  based  on  the 
electrochemical  and  thermal  modeling,  has  been  successfully 
developed  and  simulated  in  the  VTB  simulation  environment. 
Comparisons  between  the  simulation  results  and  the  experimen¬ 
tal  data,  under  different  operating  pressures  and  temperatures, 
demonstrate  the  accuracy  and  reliability  of  the  model.  A  detailed 
parameter  analysis  of  stack  working  condition  and  cell  configu¬ 
rations  of  the  SOFC  stack  is  presented  in  this  study.  Listed  below 
are  the  conclusions  derived  from  the  present  simulation  efforts: 

•  Elevated  operating  pressures  can  improve  the  cell  perfor¬ 
mance  in  terms  of  cell  voltage,  Uf  and  efficiency.  However,  the 
improvement  of  cell  performance  is  not  linear  with  respect  to 
the  increase  of  pressure  due  to  the  fact  that  the  pressure  con¬ 
tribution  to  the  Nernst  potential  is  logarithmic  in  nature. 

•  Operating  temperature  plays  a  significant  role  in  two  impor¬ 
tant  cell  performances:  Nernst  potential  and  irreversible 
losses.  Increase  in  operating  temperature  causes  a  reduction 
in  irreversible  losses  and  also  in  the  Nernst  potential.  These 
two  have  opposite  effects  on  cell  voltage  and  its  efficien¬ 
cies.  Therefore,  cell  voltage  and  efficiency  increase  initially 
and  then  with  further  increase  in  temperature  the  benefi¬ 
cial  effect  of  lower  irreversible  losses  is  completely  offset 
by  the  decrease  in  Nernst  potential,  resulting  in  a  lower 
efficiency. 

•  The  mass  flow  rates  of  air  and  fuel  affect  the  cell  performance 
due  to  the  cooling  effects  of  different  flow  rates.  However,  the 
cell  performance  can  benefit  from  the  larger  amount  of  fuel 
because  more  hydrogen  can  be  oxidized.  The  overall  influence 
of  fuel  increase  is  controlled  by  the  above  two  effects. 

•  Under  the  assumption  of  fixed  fuel  input  temperature  of  the 
SOFC,  the  mean  solid  temperature  and  hydrogen  concentra¬ 
tion  of  the  inlet  fuel  will  increase  with  the  increase  in  the 
degree  of  reforming,  which  improves  the  cell  performance  in 
terms  of  cell  voltage  and  efficiency. 

•  For  the  S/C  ratio,  higher  S/C  ratio  yields  a  lower  H2  concen¬ 
tration  and  partial  pressure,  which  has  a  negative  effect  on  the 
Nernst  potential. 

•  Increase  in  cell  diameter  will  increase  the  total  power  output 
because  of  the  larger  activation  cell  area,  however,  it  yields 
an  undesirable  increase  in  the  current  path  length  and  accord¬ 
ingly  a  higher  ohmic  loss.  Therefore,  Nernst  potential  and  the 
power  density  of  cell  both  decrease. 

•  The  length  of  the  cell  has  no  effect  on  the  cell  voltage-current 
density  relation  curve,  but  the  total  power  output  increases 
because  of  the  larger  cell  reaction  area.  However,  the  diameter 
of  the  cell  has  to  be  increased  to  counteract  the  pressure  drop 
effects  for  the  larger  cell  length,  leading  to  an  undesirable 


effect  of  increase  in  the  current  path  length  and  thereby  in  the 
cell  ohmic  resistance. 
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